Painleve III and a singular linear statistics in Hermitean 

random matrix ensembles I. 



Yang Chen 
Department of Mathematics 
Imperial College London 
180 Queen's Gates 
London SW7 2BZ UK 
Alexander Its 
Department of Mathematical Sciences 
Indiana University-Purdue University Indianapolis 
402 N. Blackford Street 
Indianapolis, IN 46202-3216 USA 

29-01-2009 



Abstract 

In this paper, we study a certain linear statistics of the unitary Laguerre ensembles, motivated 
in part by an integrable quantum field theory at finite temperature. It transpires that this is 
equivalent to the characterization of a sequence of polynomials orthogonal with respect to the 
weight 

w(x) = w(x,s) := x a e~ x e~ s / x , < x < oo, a > 0, s > 0, 

namely, the determination of the associated Hankel determinant and recurrence coefficients. Here 
w(x, s) is the Laguerre weight x a e~ x 'perturbed' by a multiplicative factor e~ s ^ x , which induces 
an infinitely strong zero at the origin. 

For polynomials orthogonal on the unit circle, a particular example where there are explicit 
formulas, the weight of which has infinitely strong zeros, was investigated by Pollazcek and Szego 
many years ago. Such weights are said to be 'singular' or irregular due to the violation of the Szego 
condition. 

In our problem, the linear statistics is a sum of the reciprocal of positive random variables 
{xj : j = 1, .., ,n}; YTj=i l / x j- 

We show that the moment generating function, or the Laplace transform of the probability 
density function of this linear statistics is expressed as the ratio of Hankel determinants and as an 
integral of the combination of a particular third Painleve function. 
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1 Introduction. 

It is a well known fact that the joint probability density of the eigenvalues, {xj : j = 1, ...,n} of 
any Hermitean matrix ensemble is [34] 

1 n 
p(x 1 ,...,x n )dx 1 ...dx n = J[ ( Xj - x k f Y\w{xi)dx h (1.1) 

Unm Kj<k<n 1=1 



where D n the normalization constant reads 



i r n 

D n [w] = - [I {x, - XkfJJwix^dxi. (1.2) 

Til lion 

+ l<j<k<n 1=1 

For the sake of concreteness the domain of integration is := [0, oo) n . 

Here w(> 0) is a weight supported on M + . Furthermore we suppose w has moments of all 
order, that is, 



( x j w(x)dx, jG{0,l,...}, (1.3) 



exist. 

We shall see that D n will play a fundamental role in this paper. 

It is also well known that the normalization constant defined above has two more alternative 
representations; the first of which as the determinant of the Hankel or moment matrix, 



AiM = det {Hj+kT jtk i Q 

n-l 



= U h r (1-4) 

3=0 

The quantity hj in the second equality of (1.4) is square of the weighted I? norm of monic 
polynomials of degree j, orthogonal with respect to w; 

/ Pj(x)P k (x)w{x)dx = hj5 jtk . (1.5) 

The random variable known as the linear statistics is a sum of a function of the random variables 
{xj : j = 1, ...,n}; 

n 

whose probability density function is determined by the standard formula, 



'(Q) = ~l [ p(xi,..,x n )S \Q-Y]f{xj) I d Xl ...dx n . 



(1.6) 
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The moment generating function denoted as Mf(s) (assuming f(x) > for all x £ R+ ) is the 
Laplace transform of Ff(Q) with respect to s and has the following form: 

poo 

M f {s) = / F f (Q)e~ sQ dQ (1.7) 
J o 

\n-l n , , „ rm-lfc..^) 

(1.8) 



det (^ +fc (s))^ =0 _ D n [w(.,s)} _ YYJJhjis) 



det (^ +fc (0));il D n [w(.,0)] n -=o MO) ' 



where 



and 



/»oo 

(s) := / x j w(x)e~ sf{x) dx, j€{0,l,..} (1.9) 
Jo 



(s) = J Pf(x)w(x)e- sf( - x) dx, (1.10) 



are moments and the square of the L 2 norm of the polynomials Pj orthogonal with respect to 
w exp(-sf), respectively. Therefore we see that the moment generating function is the ratio of 
the Hankel determinant generated by the 'perturbed' weight u>exp(— sf) to the corresponding 
quantity generated by the 'original' weight w. 

Because the moments depend on s, the coefficients of the polynomials Pj(z) also depend on 
s, however, we shall not display this dependence most of the time. 

Our monic polynomials are normalized so that: 

P n (z, s) = z n + pi(n, s) z^ 1 + ... + P n (0, s), (1.11) 

with Pq(z,s) := 1 and pi(0, s) := 0. 

We note here Heine's multiple integral representation of P n (z), 

1 f n 

p n{z) = —. V\(z - x j )p(x 1 ,...,x n )dx 1 ...dx n . (1.12) 

In section 2, a description is given for a pair of ladder operators for smooth weights. These will 
lead to a linear second order ordinary differential equation satisfied by P n (z) and two fundamental 
compatibility conditions valid for all z G C U oo. 

We denote the compatibility conditions as (Si) and (S2). 

The compatibility conditions are essentially a consequence of the recurrence relations; 

zP n (z) = P n+ i(z) + a n (s)P n (z) + /3„(a)P n _i(z), (1.13) 

together with the 'initial' conditions: Pq(z) = 1, and (3 P_i(z) = 0, and the Christoffel-Darboux 
formula (also a consequence of (1.13)). 
In the (1.13), 

a n (s) GR, n = 0, 1, .. 

and 

o j \ h n D n+ iD n -i 

Pn(s) = T = >0, n=l,2,.. 
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are the recurrence coefficients. 

An easy consequence of (1.11) and (1.13) is that 

Pi(n,s) - pi(n + l,s) = a n (s). (1.14) 

Taking a telescopic sum of (1.14) together with p 1 (0, s) = 0, implies 

n-l 

J^oy(s) = -pi(n,s). (1.15) 

j=0 

We refer the readers to [39] for basic facts about orthogonal polynomials. 
In our problem, 

f(x) := i, 0<x<oo, (1.16) 
and the unperturbed weight is given by the equation, 

w (x) = x a e~ x , a > 0. 

The compatibility conditions (Si) and (S2), and a combination of these, (S 2 ), produces a pair 
non-linear difference equations, satisfied by the auxiliary quantities a n and b n . See (2.16) and 
(2.17). 

The recurrence coefficients, cx n arid f3 n are ultimately expressed in terms of <x n and b n . See 
(2.9) and (2.14). 

The linear statistics (1.16) leads to the weight 

w(x, s) = w (x)e- s/x := x a e~ x e~ s/x , a > 0, s > 0. 

Such weights arise from a certain problem in mathematical physics: An integrable quantum field 
theory at finite temperature [31]. 

In the theory of orthogonal polynomials, the effect of infinitely strong zeros on the Hankel 
determinants, recurrence coefficients and polynomials themselves is of considerable interest. 

For orthogonal polynomials with weight w supported on [—1,1], the classical Szego theory 
gives a comprehensive account of the large n behavior of the recurrence coefficients and the poly- 
nomials (both outside [—1,1] and on (— 1, 1) ) if w is absolutely continuous and satisfies the Szego 
condition, 

" l I In w(x)\ 



rdx < 00. 



-1 \/i — x 2 

See [pp. 296-312, [39]], [23] and [35] regarding Szego's theory. 

However, there is a class of orthogonal polynomials discovered by Pollaczek and extended by 
Szego which is in some sense irregular. See [pp. 393-400, [39]] and [40] about this class of orthogonal 
polynomials. The Pollaczek-Szego weight behaves like 



exp r7CTj' c>0 ' 

near ±1, and consequently just violates the Szego condition. 
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We reproduce here some of the results of [40] to illustrate the irregularity. 
Associated with the weight, 



w(x; a, b) :- 



cosh[7T0(0)] 

where x := cos 0, < 6 < tt and 

, ... a cos + 6 , m 
W) ■= ■ n , a,b€R, a> b, 
2 suit' 

are the normalized Pollaczek polynomials {p n (x;a, b)}, 

J [p n (x; a, b)] 2 w(x; a, b)dx = 1. 

If x — > 1, then 

w(x;a,b) ~ 2 e^^ 1 " 7 ^, as -► 0, 

which shows that the weight vanishes exponentially at x = 1. An easy computation demonstrates 
the same behavior at x = — 1. 

The large n behavior of p n (x;a,b) are as follows: 

(a) p„(l;a,6) ~ n 1 / 4 e 2 ^ ^, 

(b) p n (x;a,6) ~ n K [x + (x 2 - 1)]", K = K(x) = ~^±L , x $ [-1, 1], 

zv^ — 1 

(c) p n (cos0;a,6) = A(0) cos [ra0 - 0(0) Inn + 5(0)] + e n (0), lim e n (0) = 0, < < tt. 

n— »oo 

where the n— independent functions A(0)(> 0), 5(0) are analytic in (0, it). 

This is to be contrasted with the large n behavior of the normalized Jacobi polynomials {p n (x)} 
associated with the weight 

{l-x) a {l + xf, a>-l, 0>-l, x€[-l,l], 
(a')Pn(l) ~ n^l\ 

{b') Pn (x) ~ [x+ Vx 2 -!]™, X^ [-1,1] 

(c') p„(cos0) = AxiO) cos[n0 + 5!(0)] + e n (0), 

where ^4i(0) and 5i(0) are functions of the same kind. 

The symbol ~ indicates that the ratio of the given quantities approaches a non-zero limit, 
while ±2 indicates that the limit is 1. This is the convention adopted by [40] and will not be used 
later. 

Our paper is a first step in study of the Pollaczek-Szego type orthogonal polynomials supported 
in infinite intervals. 

With reference to the Heine formula, we see that, 
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For the unperturbed or Laguerre weight, we have the explicit determination 



(-i)"p„(o.o ) = r( " + °; 1) , (Lis) 

1 [a + 1) 

since, 

„ r , xn G(n + l)G(n + a + 1) 

£>„[«;(., 0, a)] = Gja + T) ' (L19) 

where G(z) is the Barnes G— function that satisfies the functional relation G(z + 1) = T(z)G(z). 

In section 3, by taking derivative with respect to s on the orthogonality relations we obtain a 
pair of differential-difference equations or the Toda equations. Combining the Toda equations and 
the non-linear difference equations obtained in section 2 produce a particular Painleve III satisfied 
by a n (s), up to linear shift in n. 

The r— function for this Pm turns out to be intimately related to the Hankel determinant 

(poo \ n— 1 

/ x j+k x a e- x ~ s/x ) . (1.20) 
JO J j,k=0 

We also express the recurrence coefficients a n and (3 n in terms of the logarithmic derivative of 
Hankel determinant 

H n := s—lnD n (s), 
as 

and obtain a functional equation involving H n , H' n and H 1 ^. The resulting second order non- linear 
ordinary differential equation satisfied by H n , is recognized to be the Jimbo-Miwa-Okamoto a - 
form of our Pm- 

In section 4, we show, with the aid of the non-linear difference equations derived in section 
2, another functional equation involving H n , H n+ \ and H n -\- We call the resulting non-linear 
second order difference equation satisfied by H n the discrete a - form of Pm. 

In sections 5, the Riemann-Hilbert approach to orthogonal polynomials and the isomonodromy 
deformation theory of Jimbo and Miwa, are used to re-derive the Pm and thereby identify the 
auxiliary quantities, a n and b n , introduced in section 2, with the objects of the Jimbo-Miwa 
isomonodromy theory of Painleve equations. 

In section 6, we show that the Hankel determinant is the isomonodromy r— function in the 
sense of Jimbo and Miwa, and put into context of the general theory of integrable systems the 
identities derived in sections 2 and 3. 

As we are studying an example of orthogonal polynomials where the otherwise classical weight, 

w (x) = x a e~ x , x G E+, a > 

is perturbed by an infinitely strong zero, 

w(x, s) := wo(x) e~ s / x , s > 0, 

the natural questions of interest are about the large n behavior of the Hankel determinant, re- 
currence coefficients and the orthogonal polynomials. Such investigations will therefore provide 



6 



valuable insights into the asymptotic of the associated Painleve transcendant. These results will 
be published in a forthcoming paper [14]. 

We want to emphasize that in this paper and in its follow-up we do not claim the introduction of 
new concepts. Our main aim is to investigate a concrete important example of the linear statistics 
which leads to a strong zero at x = using the known techniques. In addition, taking this 
example as a "case study" we show how the apparatus which are used by the two communities 
- the orthogonal polynomial community and the integrable system community, match with each 
other. 

2 Ladder operators and non-linear difference equations. 

The pair of ladder operators has been known to various authors. See, for example, [3], [4], [7], [8], 
[9], [10], [33], [11], [18], [19], [32] and [37]. In fact, Magnus in [33] noted that such operators were 
known to Laguerre. 

Because the associated fundamental compatibility conditions and their use in the derivation of 
the Painleve transcendant [15], [16], [18], [20], [33] are perhaps less well known we summarize these 
findings in (2.1-2.5), {Si), (S'2) and (S' 2 ) in a form which we find particular easy to use. We note 
here that (Si), (S'2) and (S' 2 ) were also known to Magnus [33] and (S2) also appeared in [27]. 
See also [20]. 

For polynomials orthogonal on the unite circle the analogues ladder operators can be found 
in [25]. See [21] for the circular case applicable to bi-orthogonal polynomials. The compatibility 
condition in the circular case can be found [1] where it was used to obtain in explicit form of the 
Toeplitz determinant with the pure Fisher-Hartwig symbol and the discriminant of the associated 
orthogonal polynomials. 

The compatibility conditions can also be adapted to the situation where the weight has discon- 
tinuities. See [15] and [2]. 

The ladder operators are, 




(2.1) 



d_ 

dz 



B n {z)-y'{z))P n - 1 {z) 



= -A n _i{z)P n {z) 



(2.2) 



where 




(2.3) 



(2.4) 



v(z) 



= — \nw{z) 



(2.5) 



and the associated fundamental compatibility conditions are, 



B n+1 (z) + B n (z) = {z- a n )A n {z) - v'{z) 



{Si) 



l + (z- a n ){B n+ i{z) - B n {z)) = [5 n +iA n+ i{z) - f3 n A n -i{z), 



{S2) 



7 



valid for all z £ CU{oo}. See [12] for a recent derivation of the compatibility conditions. To arrive 
at the equations (2.3) and (2.4) we have assumed that u>(0) = w(oo) = 0. This is certainly the 
situation for our problem since we have assume that a > and s > . 

Combining suitably (Si) and (S2) gives an expression involving X^j=o^i( z )' Bn(z) and 
y'(z) from which further insight into recurrence coefficients may be gained. 

The equation (S 2 ) may be thought of as the first integral from (Si) and (S2). 

Although (S' 2 ) first appeared in [33] in a slightly different form, we present here a derivation 
of a version which we find useful in practice. 

Multiplying (S2) by A n (z) we see that the r.h.s. of the resulting equation is a first order 
difference, while the l.h.s., with (z — a n )A n (z) replaced by B n+ i(z) + B n (z) + v'(z) is a first 
order difference plus A n (z). Taking a telescopic sum, together with the initial conditions Bq(z) = 
A-i(z) = , produces the Lemma, 



Lemma 1 

n-l 
3=0 



B 2 n {z) + J(z)B n {z) + Y j A j (z) = l3 n A n (z)A n ^(z). (S' 2 [ 



If v'(z) is rational then so are A n {z) and B n (z). See (2.3) and (2.4). Furthermore, eliminating 
P n -i(z) from (2.1) and (2.2) it is easy to show that y(z) := P n (z) satisfies the second order linear 
ordinary differential equation, 

V"(z) - (v'(z) + y'(z) + ^B' n (z) - Bn^)f^ +Y,M*^J v( z ) = °- ( 2 - 6 ) 

Note that (S' 2 ) has been used to simplify the coefficient of y(z) in (2.6). 
The equation (2.6) can also be found in [37], albeit in a different form. 
For the problem at hand, f(x) = 1/x, x > 0, the weight and associated quantities are, 

w(x,s) = x a e- x - s/x , 

y(z) = z + s/z — ctlnz, v'(z) = 1 — s/z 2 — a/z 
V(z) - y'(y) 1 ( a s \ s 

z-y z\y y 2 J z 2 y' 

Using these we have the next Lemma. 

Lemma 2 The "coefficients" A n (z) and B n (z) appearing in the ladder operators are 

A n {z) = ~ z + % (2.7) 

B n {z) = + (2.8) 

z z A 

•2 

wdy, a n (0) = 0, 



:= _£_ f°5? 
K Jo y 

:= ~r~ r^^wdy, b n (0) = 0. 

hn-i Jo y 
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Proof. From the definitions of A n (z) and B n (z) and with the identities, 




obtained by integration by parts, we find (2.7) and (2.8). □ 

We see that at this stage there are 4 unknowns, a n , (3 n , a n and b n . In what follows we will 
show how (Si) and (S' 2 ) can be applied to obtain amongst other things the pair nondinear differ- 
ence equations involving a n and b n mentioned earlier. 

On equating the residues on both sides of (Si), we find, 

a n = 2n + l + a + a n (2.9) 
b n+ i + b n = s-a n a n . (2.10) 

Carrying out a similar calculation with (S' 2 ), gives, 

n-l 

(3 n = n(n + a) +b n + ^2aj (2.11) 

(3 n (a n + a n -i) = ns-(2n + a)b n (2.12) 
b\ - sb n = (3 n a n a n -i. (2.13) 

The upshot of these equations is that a n and n are entirely determined by a n and b n , where 
a n is simply a n plus 2n + 1 + a. 

Eliminating a n _i from (2.12) and (2.13) we have the next Lemma. 

Lemma 3 

(3 n a 2 n = [ns - (2n + a)b n }a n - {b 2 n - sb n ). (2.14) 

Therefore (2.14) expresses f3 n in terms of a n and b n , and importantly bypasses the finite sum in 
(2.11). Eliminating f3 n from (2.11) and (2.14), an expression can be found for ^j=o a j-> m terms 
of a n and b n . 

We state this in the next Lemma 

Lemma 4 

/ , x , , ns ~ ( 2n + a ) b n b 2 n - Sb n (n 1 -x 

2^ aj = -n[n + a) - b n H 2 ■ ( 2 -l 5 ) 

J=0 a « a n 

Note that because a n (0) = b n (0) = 0, (2.9) and (2.11) reduce to a n (0) = 2n + 1 + a, and 
/?n(0) = n(n + a) , respectively, which we recognize to be the recurrence coefficients of the monic 
Laguerre polynomials. 
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In summary, with reference to (2.13) and (2.14), we obtain two non-linear difference equations, 
satisfied by a n and b n , 

b n+ i + b n = s - (2n + 1 + a + a n )a n , ( 2 -16) 
(b 2 n - sb n )(a n + a„_i) = [ns - (2n + a)6 n ]a„a n _i, (2.17) 

to be iterated in n with the initial conditions, 

*>(') = V-s ^ a{2 ( f% y (2.18) 
b (s) = 0, (2.19) 

where K a (z) is the MacDonald functions of the second kind. 

We call (2.16) and (2.17) together with the initial conditions (2.18) and (2.19) the MacDonald's 
hierarchy. See also [4] for the treatment of a class of semi-classical weights. 

Solutions for a n and b n that are rational functions of are found for a = p + 1/2, p £ Z, 

since 

K p+1/2 (z) = K_ p _ 1/2 (z) = ^-e ^ kKp _ m2z)k - 

3 Toda evolution and Painleve III. 

By taking derivatives with respect to s on the orthogonality relation gives rise to Toda type 
equations. 
Because, 

oo 

2„ 



poo 

/ Ph 

JO 



h n = I P n wdy, 



we have, 



hence, 



d 



s ~T h n = ~s / —wdy, 
ds J y 



s-^-lnh n = -a n (3.1) 
as 

s— Infin = a n _i-a n . (3.2) 
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We also have, 



= — J P n P n - lW dy (3.3) 
= /in-i^-Pi(n) - /" ^—^-wdy (3.4) 

^ _ (3,, 

da n cia n 

(is (is 

= — b n+ i 

= fin ~ Pn+l + On- (3.6) 

The equation (3.6) follows from (2.9), (1.14) and (2.11) 
There is another identity involving X^=o a i : 

d n ^ 
j"=o 

which is an immediate consequence of a telescopic sum of the second equality of (3.6). 

We now show that the Hankel determinant D n , is up to scaling transformation the r— function 
of the Toda-equations. Let 

b n (s) := S ~< n+a )D n (s). 



We find by summing (3.1) 



d n—i 
s —\nD n = -J2aj, (3.8) 
S j=o 



since 



n— 1 

In hj = In D n . 

Applying s^ to (3.8) and keeping in mind (3.7), (2.11) and (3.8) gives 

d ( d , „ 
s — s — lnL> n 
as \ as 



/3 n - n(n + a) - a 

j=0 



n— 1 



= n - n(n + a) + s-^- In D n . 

as 



The last equation simplifies to 



In-Dn(s) = Dn+ ^ n 1 - n (n + a), 



11 



since 



Ai+1 At-1 
D 2 



In terms of D n (s) we have, 



ds 2 



In £>„(«) 



D 2 



(3.9) 



The equation (3.9) is the Toda molecule equation [38] and shows that D n (s) is the corresponding 
r— function of the Toda equations (3.2) and (3.6). 

As the Hankel determinant is now identified with the r— function (see section 6 for more on this 
issue), we may expect the emergence of a Painleve equation. In fact, a n {s) satisfies a particular 
Pin- To see this, we first investigate the evolution of a n and b n as functions of s. 

Lemma 5 For a fixed n, the auxiliary quantities a n and b n satisfy the following coupled Riccatti 
equations: 



s^- = 2b n + (2n + 1 + a + a n )a n - s 
ds 



db n 
ds 



-{b n - sb n ) + (2n + a + l)b n - ns . 



(3.10) 
(3.11) 



Proof. The equation (3.10) follows from applying to (2.9) together with the first equality of 
(3.6) and with (2.10) to replace b n+ \ by s — a n a n — b n . 

A little bit more work is required to prove (3.11). First apply s-^ to (2.11), 



d(3 n db n d ^ 
s— — = s— h s— } aj 

ds ds ds 

3=0 



db n r 
S— b n 

ds 
bl - sb n 



Sbn 



ns — (2n + a)b n — 



where the last three equalities follow from (3.7), (3.2), (2.13) and (2.14). 
The next theorem identifies particular third Painleve function. 

Theorem 1 For a fixed n G {0,1,2,..} the auxiliary quantity a n satisfies 



< = 



a n i /r> i -i i „ \ a n , a n , a 

h (2n + 1 + q)- t H — 7T H — 

s s z s z s 



1 



with the initial conditions 



□ 



(3.12) 



o n (0) = 0, a' n (0) = -, a>0. 



a 



(3.13) 



// a n (s) := — q(s), then q(s) is Pnp(— 4(2n + 1 + a), — 4a, 4, — 4), following the convention of 
[36]. 
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Proof. Eliminate b n from (3.10) and (3.11) gives (3.12). The initial conditions follows from a 

straightforward computation. □ 

Remark I. If n = 0, then bo = and (3.10) is solved by 



K a+l {2^s) 



which is (2.18). We observe that the above also solve (3.12) for n = 0. 
Remark II. An alternative form, obtained by, 



reads, 



VY„ — 



K 

s 



a n (s) 



X n (sY 



2n + 1 + a X% 



X n s s 2 s s 2 X n 



(3.14) 



which is a Pnr{— 4a, — 4(2n + 1 + a), 4, —4). If the derivatives in (3.14) were neglected, then X„ 
solves the quartic, 



X 4 - aX 3 - (2n + 1 + a)sX - s 2 



0. 



(3.15) 



We may interpret an appropriate solution (3.15) as the geometric mean of the end points of the 
support of a single interval equilibrium density. This appears in a potential theoretic minimization 
problem, the detail of which is in a forthcoming paper [14]. We note that another Pm associated 
with the Toeplitz determinant 

det[7 j _ fe+J ,(v / t)]o<j ,k<n—li 

was found in [21] (see also [22]). Here I r (z) is the modified Bessel's function of the first kind. 
We note that this Toeplitz determinant can be thought of as a Toeplitz-analog of the Hankel 
determinant we are studying in this paper. Also, the above Toeplitz determinant with v = 
appeared [41] in connection with a certain ensemble of n x n unitary matrices and Ulam's problem 
in combinatorics. 

In the next theorem, we display two alternative integral representation of D n , in terms of a n 
and X n . 



Theorem 2 



hi 



Dn{s) 

D n (Q) 



t 

(In 



- f 

JO 



( a\ a 2 1 / 
-°»("+2)-f + l( 

!-i<*.-°> j -(«+!) 



/2 



X n AX 2 AX 2 



dt 

T 



dt 

T 



(3.16) 
(3.17) 



Proof: From (3.8) and (2.15) we see that the logarithmic derivative of D n {s) is expressed in terms 
of a n and b n . If we use (3.10) to eliminate b n in favor of a n and a! n from the resulting equation, 
then (3.16) follows after some simplification. The equation (3.17) follows from the substitution 
a n (s) =: s/X n (s). Note that D n (0) is given by (1.19). □ 
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Put 

H n :=s^-lnD n . (3.18) 
as 

In section 6 we will show that the Hankel Determinant D n (s) can be identified with the Jimbo- 
Miwa r -function corresponding to the solution a n (s) of the Painleve III ' equation (3.12). Namely, 
we will show that 

- n ( n + Q ) 

D n (s) = const r(s) eas 2 . (3.19) 
In its turn, this relation yields the following formula for the quantity H n , 

, . s n(n + a) . 
Hn = ct(s) + - + 1 - ; , (3.20) 

where a(s) = s^lnr is the Jimbo-Miwa-Okamoto a -function corresponding to the equation 
Pi a'. Therefore, H n , H' n and H'^ should satisfy a functional equation 

f(H n ,H' n ,HZ,n,8,) = 0, 

known as the Jimbo-Miwa-Okamoto a - form of our Pup. 

With H n defined above, it is easy to see from (3.7), (3.8) and (2.11) that 

b n = sH' n (3.21) 
Pn = n(n + a) + sH' n - H n . (3.22) 

In the next theorem we state the non-linear second order ordinary differential equation satisfied by 



Theorem 3 // 



H n := s^- In D n (s), (3.23) 



then 

{sH'^) 2 = [n- (2n + a)H' n f - 4[n(n + a) + sH' n - H n ]H' n (H' n - 1). (3.24) 
Proof. First we re-write (2.14) and (3.11) as 

(3 n a n + 6 " ~ Sbn = ns-(2n + a)b n (3.25) 
— {b 2 n -sb n ) = sb' n - b n + ns - (2n + a)b n (3.26) 

respectively. Eliminate a n from (3.25) and (3.26) produces 

« - b n ) 2 = [ns - (2n + a)b n f - A(3 n {b 2 n - sb n ). (3.27) 
The equation (3.24) follows by substituting b n and (3 n from (3.21) and (3.22) into (3.27). □ 
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Hence the recurrence coefficients a n and (3 n , of the orthogonal polynomials associated with 
our weight, 

w(x,s) = x a e~ x ~ s/x , x£[0,oo), a > 0, s > 0, 
are expressed in terms of H n , H' n and H", as follows: 

2s(H' ) 2 - ?H' 

a n = 2n + l + a + „ , \ , (3.28) 

n = n{n + a) + sH' n -H n , (3.29) 
and that H n itself satisfies (3.24). 

Remark III. With the identification (3.20) of the function H n as the a -function (up to a linear 

shift), equation (3.24) coincides, up to the changing the independent variable s to the variable 
t = ^fs, with equation (C.29) of [30] 1 . 



4 Discrete a— form 

We may anticipate due to the recurrence relations and non-linear difference equations, (2.16) and 
(2.17) that, for a fixed t, H n , and H n ±\ would satisfy a discrete analog of (3.24). It turns out 
that in this instance a n has a simpler expression in term of H n ±±. We shall obtain an expression 
for b n in terms of H n and H n ±\. From 

n-l 

H n = - y^aj, 
j=o 

we have 

a n = H n — H n+ i, (4.1) 

and 

«n + «n-i = H n -\ — H n+ i =: 5 2 H n . 

Multiply the above by (3 n we find, 

/3 n (a n + a„_i) = l3 n 5 2 H n (4.2) 
= ns- (2n + a)b n , (4.3) 

where the last equation follows from (2.12). Now (2.11) becomes 

P n = n(n + a) +b n - H n . 



1 When comparing equation (3.24) with equation (C.29) of [30] one has to also take into account that, if 
we denote the a -function of [30], for which equation (C.29) is written, as <jjm then the exact relation with 
our a function is given by the equation, <JjAi(t) = 2a(t 2 ). 
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Substituting the above into (4.2) produces a linear equation in b n whose solution is 

_ ns + 5 2 H n [H n - n(n + a)} 
n ~ 2n + a + 5 2 H n ' 



(4.4) 



Hence the auxiliary quantities a n and b n , and (3 n are now expressed in terms of H n , and H n ±\. 
Substituting these into (2.13) give rise to the discrete a— form stated in the next Theorem. 

Theorem 4 // 

H n := s— lnD„, 
as 

then 

{[H n - n(n + a)]5 2 H n + ns}{[H n - n(n + a) - s]5 2 H n - (n + a)s} 
= {2n + a + 5 2 H n ){ns + (2n + a)[n(n + a) - H n ]}(H n - H n+1 )(H n ^ - H n ). (4.5) 

We also have the discrete analog of (3.28) and (3.29), 

a n = 2n + 1 + a + H n - H n+1 (4.6) 

ns-n(n + a)5 2 H n -(2n + a)H n 
Pn = »(" + «) + 2» + a + M w • (47) 

Now the obvious equalities, (4.6) = (3.28) and (4.7) = (3.29), imply two further differential- 
difference equations which must satisfy, 

2s(H' n ) 2 - sH' n 
Hn ~ Hn+1 ~ sH» + n- (2n + a)H' n (48) 
ns-n(n + a)5 2 H n -(2n + a)H n , 

2n + a + 5 2 H n = ~ (49) 



Remark IV. From the point of view of general Jimbo-Miwa isomonodromy theory of Painleve 
equations, which we will be discussing in sections 5 and 6, equation (4.5) should be related to the 
Backlund - Schlesinger transformations of the r - function. However, we failed to identified equation 
(4.5) with any of the difference equations for the r - function discussed in [30] and describing the 
possible Schlesinger transformations. Since it is written for the logarithmic derivative of the r 
- function, and not for the r - function itself as in [30], equation (4.5) might be in fact of a 
different nature than the ones considered in [30]. We also want to mention that, equation (4.5) 
is an integrable discrete equation - its Lax pair is formed by the first and the third equations of 
the triple (5.17) of section 5, and the equation itself, as it has already been noticed, represents a 
Backlund - Schlesinger transformation of the third Painleve equation. Therefore, we expect this 
equation to be equivalent to one of the known discrete Painleve equations, which we have not 
yet identify. Apparently this identification is not quite straightforward. One of the referees has 
suggested that (4.5) could be a composition of the basic Schlesinger transformation T\ and Ti 
found in [22] (see Proposition (4.6) of [22]). 
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5 An alternative derivation of the Painleve III equa- 
tion. 

In this section we present an alternative derivation of the third Painleve equation (3.12) for the 
quantity a n (s). This derivation is based on the Riemann-Hilbert point of view [18], [19] on or- 
thogonal polynomials and makes use of the general Jimbo-Miwa-Ueno theory of isomonodromy 
deformations. This in turn allows us to place some of the key identities of the preceding sections 
into a general framework of integrable systems. 

The Riemann-Hilbert problem for the orthogonal polynomials at hand is the following 

• Y : C \ R+ -> C 2x2 is analytic. 

/ ' ~| ft — x — s f x\ 

• y+(x) = Y-(x) ( X e 1 J for x € K+ \ {0} , with M+ oriented from left to right. 



1 



. V!:! (1 • Oil :))(^' . "„ I as : - x . 
• Y(z) = 0(1) as z -> 0. 



Here Y±(z) denote the non-tangential limiting values of Y(z) on R + taken (in the usual point- 
wise sense) from the ± - side. The Riemann-Hilbert problem has the unique solution expressed in 
terms of the orthogonal polynomials P n (z) , 



( 



Y(z) 



Pn(z) 



^ Jr. 



K- 



P n {x)x a e- X - S ' x ^ \ 

x — z 
P n ^i(x)x a e~ x ~ s l x 



dx 



x — z 



(5.1) 



7 



We also note that an immediate consequence of the uni-modularity of the jump matrix of the above 
Riemann-Hilbert problem is the identity 2 , 



detY(z) = 1. 



(5.2) 



Equation (5.1) implies, in particular, that the asymptotic behavior of the function Y(z) at 
z = oo and z = can be specified as the following full asymptotic series, 



Y, 



OO, 



k=l 



and 



Y{z)~Q (l + J2Y k z k J, z^O, 



(5.3) 



(5.4) 



k=i 



2 A more conventional derivation of this identity is based on the use of the basic three-term recurrence 
equation (1.13) (see e.g. [13]) 
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where 0-3 denotes, as usual, the third Pauli matrix 

'l 



0-3 








Moreover, the matrix coefficients Y±k and Q of these series are the smooth functions of n and 
s (and of course of a), and they all can be easily expressed in terms of the fundamental objects 
associated with the orthogonal polynomials P n (z) , i.e. in terms of the functions p^.(tx, s), /^(s), 
and the negative moments of the polynomials P n (z) . Indeed, by a straightforward calculation 
we have from (5.1) the following expressions for the coefficient F_i and for the matrix multiplier 
Q = Y(0) . 

Pi(") ~ti , 

(5.5) 
(5.6) 

\-q n 1-pnQnJ 

where 

P n (0) f P n {x)x a e- X - S ' x A 2m P n _i(0) 

Pn = ^— / ds, g„ = , 5.7) 

2m J R+ x hn-i P n (0) 

and we have taken into account the determinant identity (5.2). 

We are now going to write down the triple of the differential and difference equations for the 
function Y n (z,s) following the standard procedure of the theory of integrable systems (see [17], 
[29], [30]; see also [18], [26], [13], and [28]). 

Put 

V(z) = V{z,n,s,a) := y(z)e - H*+i ) a3 z% a3 , (5.8) 

where the branch of the function z? we define by the condition, — tt < argz < it . The Riemann- 
Hilbert relations in terms of the function ^f(z) reads as follows, 

(i) * : C \ R -> C 2x2 is analytic, 
(ii) *+(x) = *_(z) ^ for x G R and x > . 

(iii) *+(x) = <&-(x)e niaa3 for x G R and x < . 

(iv) ~ (j + Er=i %r) z( n+ ^) CT3 e-f CT3 asz^co. 

(v) ~Q(/ + Er=i^^') zfa3e ~^ a3 as *->°> 

where the real line R is oriented as usual from left to right. The coefficients *&±k of the asymptotic 
series are easily evaluated via combinations of the coefficients Y±k and the coefficients of the 
expansions of the exponential function e~5i and e~2 near z = 00 and z = 0, respectively. In 
particular, we have that, 

>w-§ -fe , 

*-i = • (5-9) 

-Pi(») + * 
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The important feature of the -RH problem is that the jump matrices of the jump relations (ii) 
and (iii) do not depend on z, s, and n . Therefore, by standard arguments based on the Liouville 
theorem (cf. [17], [18]), we conclude that the logarithmic derivatives, 

A t z y-!Wf±y-U z \ b(z) :=^^-^- 1 (z), and U(z) := tffz, n + l)^" 1 ^, n), (5.10) 

OZ OS 

are rational functions of z. Using the asymptotic expansions (iv) and (v), we can evaluate the 
respective principal parts at the poles at the points z = and z = oo and arrive, taking into 
account (5.9) and (5.6), at the following explicit formulae for the function A(z) , B(z) , and U(z) . 



and 



where 3 



A(z) = --a 3 + - 



B{z) = - 



A, 



sz 




U(z) = 

A 1 = I[ < 7 3 ,tf_ 1 ] + (n + T)a 3 



hn-l 



2m 

—n — i 



(5.11) 
(5.12) 

(5.13) 



(5.14) 



and 



A 2 = fazQ' 1 = | 



1 - 2p n q n -2p n 
2q n (p n q n -l) 2p n q n -\ 



1 0\ II 0^ 
C/ = *-i(n + l) | - I | tf-i(n) 

,0 0/ 



,0 0, 



'pi(n + l)-pi(n) ti 



2tH 





f—a n 


h n \ 

2ni \ 


H 


V h n 





(5.15) 



(5.16) 



According to the standard methodology (cf. [29], [30], [18]), relations (5.10) should be now 
re-interprited as a system of linear differential-difference equations, 



' *T = A(zMz) 
9 -W = B(zMz) 
#(z,n + 1) = U(z)W(z, n), 



(5.17) 



3 Notation [Mi, M2] means the usual commutator of two matrices, [Mi, M2] = M1M2 — M2M1 . 
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which we call the Lax triple. 

The compatibility conditions of this system, i.e. the equations, 

^-^ = [^MW] = (5.18) 

^1 =B{z,n+l)U{z)-U(z)B(z,n), n + l)) s = tf a (z, n + 1)) (5.19) 



and 



dU(z) 



_ =A(z,n+l)U(z)-U(z)A(z,n), n + l)) z = V z {z, n + 1)) (5.20) 

yield the Painleve type (equ. (5.18)), the Toda type (equ. (5.19)) and the discrete Painleve or 
Freud type (equ. (5.20)) equations, respectively, for a proper combination of functions h n , p n and 
q n . Moreover, using the Jimbo-Miwa list of the Lax pairs for Painleve equations [30] and noticing 
that the "master equation", i.e. the first equation of system (5.17) is a 2x2 system with two 
irregular singular points of Poincare rank 1, one concludes that the relevant Painleve equation is, 
in fact, the third Painleve equation. In order to make a precise statement, i.e. to point out the 
exact combination of functions h n , p n and q n which makes the solution of Painleve III equation, 
we only need to perform a simple scaling transformation of system (5.17) which would bring it to 
the normal form of [30]. To this end, we introduce the new independent variables, 



\:=s' 1/2 z, and t := y/s, (5.21) 
and pass from the function s) to the function 3>(A,i) defined by the equation, 

$(A, t) := H n+ ? >^(tX, t 2 ). (5.22) 

We notice that in terms of the function <b(\,t) the asymptotic relations (iv) and (v) transform 
into the relations, 



oo 



*(A)~ ^ + E^?J A(n+f)fT3e "^ 3 ' A -°°> ( 5 - 23 ) 
and 

$(A) ~R (^ I + Yj^k^ Af^e-^ 3 , A^O, (5.24) 
with the new coefficients connected to the old ones by the equations, 

Simultaneously, the first two equations of system (5.17) transform into the Jimbo-Miwa Lax pair 
for the third Painleve equation, 



^ = ("1-3 + ^ + ¥) - A JM (XMX) 

(5.25) 



a$(A) 

dt 



(-|a 3 + Bq + $(A) = Bjm(XMX). 
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Here, the matrix coefficients A_i , A_ 2 , Bq , and B_i are given by the equations, 

u 



A-i = H n+ f ) CT3 Ait( n+ f > :i = [ 2 ) , (5.26) 

v 



a_ 2 = V( n+ t)°M 2 t( n+ f )°* 



C+3 -< 
C+t /■ 

w ^2, 



(5.27) 



where 



^o = -^- 1 - i (n + -).3 = -^ J, (5.28) 



i i /c + 1 -< 

5-i = ~ 7 A-2 = -t I , (5.29) 

\ w ^> 2, 



oo = - a - 2n, (5.30) 



and the new scalar functional parameters u , v , £ , and w are defined in terms of the original 
functions h n , h n -i , p ra , and q n via the formulae, 

2m h n -i 

C = -tp n q n , w = -—r 2n - a . (5.32) 
q n 

With the u — w notations, the Lax pair (5.25) matches, up to the replacement t —t and the 
use of the letter z instead of the letter £ , the Lax pair presented on page 439 of [30] , and hence 
we can use the general results of [30]. 

Theorem 5 ([30]) Consider the over- determined linear system (5.25) with the matrix coefficients 
defined by the equations 





<--¥ u ) 




H 








K v e -tJ 





C + I ~w( 



C+t /■ _ t 



(5.33) 



v££ —C 



w s 2, 



i (o u\ 1 k + k -< 

= T , B_i = -- | , (5.34) 

i.e., by the right hand sides of the last equalities in the formulae (5.26) - (5.29) not necessarily 
assuming any specific, "orthogonal polynomial" choice of the parameters 9^, u, v, C and w . Then 
the following is true. 
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1. The compatibility condition of system (5.25), i.e. the matrix equation, 

9Ajm 9Bj M 



dt d\ 

is equivalent to the following set of scalar equations, 



= [Bjm,Aj M }, (5.35) 



with the quantity, 



du 

t— = Ooou + 2t(w, (5.36) 

t^ = -9coV + -(t + t), (5.37) 
at w 

4 = 2( » W + C+ MC±^), (5.38) 
at w 

d\nw 2u „ 

t—j— = 2WV-0OO, 5.39 

at w 

(b :__^ (2f + ()+ *iS±a_^ TO , (5.40) 



being the first integral of system (5.36) - (5.39). 

2. The quantity 6q is very similar to the parameter 6^ . Indeed, they describe the formal 
monodromy at the relevant irregular points; 6^ - at A = oo and 9q - at A = 0. This 
means that they appear as the branching exponents in the following formal matrix solutions 4 
of the Lax pair (5.25) at the points A = oo and A = . 

/ °° ^ (00) \ 



and 



3. The function 



*trnjn = Ri0) [l + Y,*k )xk ) xU °-* n - ^ 



k=i 



a 

V ■= -y-, (5.43) 
(,w 



satisfies the third Painleve equation, 

In the standard notations [24], this is the Painleve III equation, Pni(— 4#o, — 4(1— 9oo), 4, —4) . 



4 All the coefficients < 5>k°°' ) of series (5.41) are uniquely defined as rational functions of u, v, C,w, and t 
via simple recurrence relations [29]. The matrix factor is defined, up to a right diagonal multiplier, 

R(0) r(0) A ) by thc equations, fi? (0) cr 3 (R^y 1 = A_ 2 , det = 1 . The coefficients $j; 0) of series 

(5.42) are defined as rational functions of u,v,Q,w, and t up to the conjugation — > A _1 <I>j^A. We 
will discuss this issue in more detail in the next section. 
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Comparing the general formal expansions (5.41), (5.42) with the asymptotic series (5.23), (5.24), 
we conclude that in our case the parameters 6oo and 6o assume the values (see also (5.30)), 

#oo = —a — 2n, and 9$ = «• (5.45) 

Simultaneously, from (5.31)-(5.32) we see that in our case, 

= (546) 
2vrztp„ V ' 

In other words, relation (5.46) provides us with the combination of the functions h n ,p n ,q n which 
satisfies the third Painleve equations and which we have been looking for, while relations (5.45) 
specifies the parameters of the Painleve equation. We conclude then that the function y defined 
in (5.46) satisfies the following Painleve III equation, 

§-K*H£-^ + «* +i+ °>) + v-j. (5 - 47) 

What is still left for us is to establish the connection between the function y{t) defined in (5.46) 
and the function a n (s) defined in (2.7). This is easy; indeed, we have that 

J -JL w dx = J P(x) (x n ' 1 + p 1 (n)x n - 2 + ... + ^p-\wdx = P n (0) J -j^wdx, 

and hence (see (5.7)), 

2iri s ' 



Pn = 7T- — > ( 5 - 48 ) 



which in turn implies that 

V = - = ^. (5-49) 
It is an elementary exercise to check that the substitution, 

y = — , t = y/s, 

transforms equation (5.47) into equation (3.12) for the quantity a n (s) . 

We conclude this section by revealing the connections of some of the key identities established 
in the previous sections with the general constructions of the isomonodromy theory of Painleve 
equations discussed above. 

We first notice that, similar to the derivation of equation (5.48), the following 

q n = . (5.50) 

Therefore all the matrix coefficients of the Lax triple (5.17) can be expressed in terms of the 
functions a n ,h n ,a n , and b n . Actually, we need only to re-write the coefficient A2 from (5.15), 

A 2 = - I . (5.51) 



2 



4-7U b n ( bn _ -A 2b n _ 2 
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Substituting this representation together with the similar formulae for A\ and Uq (see (5.14) 
and (5.16)) into the compatibility equation (5.20) yields the set of the scalar difference equations 
(2.9) - (2.13) of section 2. Similar operation with the compatibility equation (5.19) results in 
the differential-difference equations (3.1) - (3.6) of section 3. The scalar form of compatibility 
equation (5.18) we have already discussed in detail in this section. As we have seen, this equa- 
tion is equivalent to the set of scalar equations (5.36) - (5.39) which, in terms of the functions 
h n (s), /i ra _i(s), a n (s), b n (s) , transforms to the system, 

S ~Js = ~ hnan > (5.52) 

^ = ^H-y, (5.53) 
ds h n a n 

= jj n _ hl( s _ 5 n ) _ -^-a n (5.54) 
ds a n h n -\ 

S ~d~s~ = 2 ^ n + + 1 + 01 + an ^ an ~ s ' (5.55) 
which is equivalent to the system of equations derived in Lemma 5 of section 3. 

We want to highlight the theoretical meaning of the important identity (2.14). It is, in fact, the 
formal monodromy identity (5.40) written in terms of the functions h n (s), h n ~\(s), a n (s), b n (s) . 

It is also worth mentioning, that the ladder equations (2.1) and (2.2) are just the first column 
of the master, z - equation of the Lax triple (5.17). The first column of the second equation of the 
triple (5.17) yields the ladder operators in s , 

d 



(zs— - b n ^P n (z) = -f3 n a n P n _i(z), (5.56) 



ZS ~(L ~ bn ~ 1 ~ (Xn - ian - 1 + za n-\j P n -\{z) = a n - 1 P n (z). (5.57) 

When deriving these equations we made use of the relations (2.10) and (2.13). The equations (5.56) 
and (5.57) of course can be derived using the orthogonality relations, bypassing the isomonodromy 
theory. 

Finally, the third equation in (5.17) is equivalent to the basic recurrence relations (1.13) for 
the polynomials P n . It should be pointed out, that for the general solutions of the Painleve III 
equations, only the first two equations, which are equivalent under scaling to the Lax pair (5.25), 
take place. From the point of view of the general isomonodromy theory of Painleve equations, the 
third equation of the triple (5.17) describes the Backlund-Schlesinger transformation of the third 
Painleve equation (see [29]; see also Chapter 6 of [28]). 

6 The Hankel determinant as the isomonodromy r - 
function. 

Let us remind the reader the Jimbo-Miwa definition of the r -function corresponding to the third 
Painleve equation (5.44). 
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Consider the formal series (5.41) and (5.42) of theorem 5. In the footnote to this theorem we 
have already mentioned that all the coefficients of these series can be evaluated as rational functions 
of the parameters u, v, £, w, and t . For example, the first coefficients in each of he series are given 
by the relations (see p. 440 of [30]), 



+2 



^ = 7 J, (6.1) 



-uv + t( + \ 



and 

i 



i (uv — tC — £ —u 
<>=7 , I • (6.2) 



t 



v -uv + tC + V/ 



2 

where the parameters u and v are related to the basic parameters u and v through the equation, 

*°yV (!_♦)■ 

which, in particular, means the identity, 

n2 _ n2 

uv = uv + ^2£ — Zo. ( 6 .4) 
4 

Denote Y'oo(A) and ^o(A) the series in the brackets of formulae (5.41) and (5.42), i.e. 

/ oo $(oo)\ 

no(A)=^ + X;^E-J. (6-5) 



and 



oo 



Y (X)=^I + ^f ) X k j. (6.6) 

The Jimbo-Miwa-Ueno isomonodromy r -function [29], in the case of the Lax pair (5.25) is defined 
by the formula, 

BY cfY 
dlnr = -TraceRes A=0 y - 1 (A)^(A)dT (A) - Trace Res A=oo y^ 1 (A)^(A)dT 00 (A), (6.7) 

where ^ 

dT (X) = -— a 3 dt, and dT^X) = --a 3 dt. 

Substituting (6.5) and (6.6) into (6.7), we arrive at the equation, 

dlnr = ^ Trace ($f ] (J 3 + ^a^j dt, 

which, taking into account (6.1), (6.2), and (6.4), implies that, 

dkn-(t) = H IH (u(t),v(t),((t);t)dt, (6.8) 



25 



where 



H in (u, v, c:/) | Uuv - 2t( I- ^_^c j , (; f)) 



Let us now turn to the Hankel determnant D n (s) and consider the quantity 

H n := s—lnD n , 

as 

which has played a central role in sections 3 and 4. From (3.8) and (2.15) we have that 

ns — (2n + a)b n b 2 — s6„ 
H n = n(n + a) + b n ^ — + " 2 - ■ (6.10) 

At the same time, in our case (see (5.45)), 

b n 

6 = a, 8 OQ = -a-2n, uv = -[3 n , and ( = — — , 

where for the last two formulae we have used (5.31), (5.32), (5.48) and (5.50). Substituting these 
formulae into (6.9) we obtain that for our special solution of system (5.36) - (5.39) the function 
Hju assumes the form, 

Hin = - (-2f3 n + 2b n -s + n(n + a)) . (6.11) 

Recalling now identity (2.14) (which we remind the reader is the formal monodromy identity (5.40) 
in disguise) we arrive, after some simple algebra, at the relation, 

H HI = \ Hn -t- n ±±^. (6.12) 

This equation, in its turn, implies that 

d t d t s n(n + a) 

s— Inr = -— inr = -Httt = H n 

ds 2dt 2 2 2 

d s n(n + a) 

= S Ts lnDn -2-^2—' (6 - 13) 

and hence we obtain the following relation between the Hankel determinant D n and the r -function 

of the third Painleve equation (compare with [26] where similar formula is derived for a class of 
Toeplitz determinants), 

s n(n-f-CK) 

D n (s) = constr(s)e2s 2 . (6.14) 

Remark V. The function Hjjj(u,v, (;t) actually depends on ( and the product uv . The latter, 
can be expressed, with the help of the formal monodromy relation (5.40) in terms of £ , t , and y . 
Indeed, we have that 

2uv = (6 + e^ty + 2Cy#oc + 2y 2 C(C + *), 

and expression (6.9) can be transformed to the following equation defined Hju as a function n y , 
£ and t , 

Hm(y, (;t) = - t (VC 2 + 2(ty 2 + O^y - t)( + (0 O + O^ty - t 2 - ^-^) • (6.15) 
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This is the canonical (see again p. 440 of [30]) representation of the logarithmic derivative of the r 
- function for Painleve III equation (5.44). The remarkable fact of the general Jimbo-Miwa-Ueno 
theory is that the function Hjjj(y,(',t) is the Hamiltonian of the third Painleve equation. 

Remark VI. In sections 2 and 3, an important role has been played by equation (2.15) which 
transforms a non-local object - the sum of aj from j = to j = n — 1 , to a local expression, which 
involves only a n and b n . We can see now an intrinsic reason for that. Indeed, on one hand, the 
sum mentioned is, by its very nature, the logarithmic derivative of the Hankel determinant. On 
the other hand, the latter is a r - function and hence its logarithmic derivative must admit a local 
representation in view of the general formula (6.7). 

Remark VII. The isomondromy context for orthogonal polynomials and, in particular, the inter- 
pretation of the Hankel determinants as isomonodromy r - functions have been well understood for 
some time, since the early nineties works [18] and [33]. For the most general semiclassical weights 
this fact was established in the recent paper [5] (see also [6]) 
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